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ABSTRACT 

A transient  ballistics  and  combustion  model  is  derived  to  represent  the 
closed  vessel  experiment  that  is  widely  used  to  characterize  propellants.  The 
model  incorporates  the  nitramine  combustion  mechanisms  which  are  contained  within 
the  steady-state  model  developed  as  a part  of  last  year's  (FY  1976)  research.  A 
computer  program  is  developed  to  solve  the  time-dependent  equations,  and  is 
applied  to  explain  aspects  of  closed  vessel  behavior.  It  is  found  that  the  rate 
of  pressurization  in  the  closed  vessel  is  insufficient  at  pressures  of  interest 
to  augment  the  burning  rate  by  time-dependent  processes.  In  the  case  of  nitramine 
propellants,  however,  the  cratering  of  the  burning  surface  associated  with  combustion 
above  break-point  pressures  augments  the  effective  burning  rate  as  deduced  from 
the  closed  vessel  experiment.  Low  pressure  combustion  is  significantly  affected 
by  the  ignition  process  and,  in  the  case  of  nitramine  propellants,  by  the  developing 
and  changing  surface  structure.  Thus,  burning  rates  deduced  from  the  closed  vessel 
experiment  may  or  may  not  agree  with  those  measured  in  the  equilibrium  strand 
burner.  Application  of  the  model  to  closed  vessel  test  cases  shows  good  agreement 
between  theory  and  experiment. 

Series  of  T-burner  experiments  are  performed  to  compare  the  combustion  instability 
characteristics  of  nitramine  (HMX)  - containing  propellants  and  ammonium  perchlorate 
(AP) propellants.  It  is  found  that  the  inclusion  of  HMX  consistently  renders  the 
propellant  more  stable.  Although  ash  produced  by  more  fuel-rich  propellants  could 
have  provided  mechanical  suppresssion,  results  from  clean-burning  propellants  permit 
the  conclusion  that  HMX  reduces  the  acoustic  driving. 

Additional  strand  burning  rate  data  are  presented  which  provide  further  confirma- 
tion of  the  results  of  the  steady-state  modeling  effect  of  FY  1976. 
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SECTION  1 
OBJECTIVE 


The  objective  of  this  research  program  Is  to  study  two  aspects  of  the  transient 
combustion  properties  of  nitramine  smokeless  propellants.  One  is  to  determine  the 
role  of  the  nitramine  ingredient  in  the  combustion  driving  of  acoustic  instability 
in  solid  rocket  motors.  The  second  is  to  determine  the  effect  of  a transient 
pressure  environment  on  the  combustion  of  nitramine  propellants  in  armament  or 
closed  vessel  apparatus.  The  steady-state  combustion  characteristics  of  nitramine 
propellants  have  been  described  as  a result  of  work  performed  under  AFOSR  Support 
Agreements  AF0SR-ISSA-75-0005  and  AFOSR-ISSA- 76-0006  (Ref.  1).  The  present  research 
begins  to  investigate  unsteady-state  or  transient  combustion  characteristics. 
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SECTION  2 
INTRODUCTION 

Nltramlnes  are  of  Interest  for  gun  and  rocket  propellant  applications  because 
they  are  an  energetic  source  of  smokeless  combustion  products  (Refs.  2-4).  Recent 
studies  have  been  devoted  to  understanding  and  Improving  the  steady-sitate  combustion 
characteristics  of  nltramlne  propellants  (Refs.  1,  5-10).  Rocket  motors,  are, 
generally,  principally  concerned  with  steady-state  performance.  Although  It  is 
recognized  that  gun  propellants  perform  in  a highly  transient  pressure  environment, 
there  Is  basis  for  the  Inference  that  Improvements  In  steady-state  combustion  will 
be  reflected  as  Improvements  in  the  transient  ballistics  performance  (Ref.  11). 

Anomalous  combustion  behavior  can  be  brought  about  by  unsteady  processes  occur- 
ring In  combustion  chambers.  In  rocket  motors,  the  most  prevalent  behavior  is  com- 
bustion instability.  Combustion  instability  Is  of  particular  concern  in  smokeless 
propellants  because  the  absence  of  aluminum  powder  deprives  the  system  of  significant 
stabilization  (Ref.  12).  Virtually  no  work  has  been  done  to  describe  the  influence  of 
nitramlnes  on  combustion  instability  properties.  In  projectiles,  the  extreme  trans- 
ient pressure  environment  can  alter  the  combustion  process  from  the  steady-state  as 
a result  of  time  lags  in  the  combustion  response.  Although  this  problem  has  been 
studied  extensively  with  respect  to  rocket  motor  Ignition  and  controlled  termina- 
tion (e.g..  Refs.  13-17),  it  has  not  been  applied  to  gun  propellants  or  to  the  more 
extreme  projectile  environment.  Burning  rates  derived  from  closed-vessel  (transient 
pressure)  testing  differ  from  those  measured  at  constant  pressure,  which  raises  a ' \ 

question  about  the  predictability  of  ballistics  performance  (Ref.  18).  Further,  the 
transient  pressure  environment  may  (Ref.  11)  or  may  not  (Ref.  19)  aggravate  any 
steady-state  combustion  anomaly  to  Impair  the  quality  of  the  Impetus  delivered. 

It  Is  considered  that  understanding  the  transient  combustion  characteristics  of 
nltramlne  smokeless  propellants  will  provide  additional  direction  for  future  develop- 
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mental  efforts.  The  mechanistic  understanding  derived  from  the  steady-state  research 
provides  an  essential  background,  and  an  understanding  of  the  transient  phenomena 
will  afford  Improved  tailoring  and  predictive  capabilities. 
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SECTION  3 

TRANSIENT  COMBUSTION  MODEL  DEVELOPMENT 

! 3.1  REVIEW  OF  CLOSED  VESSEL  BURNING  RATE  DATA 

A review  was  conducted  of  available  data  wherein  direct  comparisons  were  made 
between  burning  rates  as  measured  in  a closed  vessel  (transient  pressure)  and  as 
measured  in  a strand  burner  (equilibrium  pressure)  (Refs.  20-24).  Some  of  the  closed 
vessel  data  were  accompanied  by  the  associated  test  pressure-time,  or  pressure-dp/dt 
(P-dot)  histories.  As  will  be  discussed  below,  it  is  important  to  know  the  P-dot 
values  for  transient  analysis.  Some  of  the  data  were  for  nitramine  propellants; 
some  were  for  homogeneous  (unfilled  single-base  or  double-base)  active  binder  pro- 
pellants. In  general,  there  appears  to  be  very  good  agreement  between  the  closed 
vessel  and  strand  data  for  the  homogeneous  propellants.  Such  was  the  conclusion  of 
Grollman  (Ref.  23)  in  his  study.  Figures  contained  in  Ref.  (24)  show  a virtual 
identity  between  closed  vessel  and  strand  data.  Lenschitz'  data  for  homogeneous 
propellants  (Ref.  20),  some  of  which  are  reproduced  here  as  Fig.  1,  also  show  good 

. j 

agreement  between  the  two  methods.  On  the  other  hand,  there  is  a significant  and 
consistent  difference  between  closed  vessel  and  strand  burning  rates  in  the  case  of  { 

propellants  containing  HMX  (Refs.  20-22).  Examples  from  Ref.  (20)  are  presented 
here  as  Fig.  2.  At  low  pressure,  the  closed  vessel  rate  tends  to  lag  the  strand 
rate;  at  high  pressure,  the  closed  vessel  rate  tends  to  lead  the  strand  rate.  As 
a result,  the  closed  vessel  data  exhibit  a more  uniform  and  higher  apparent  pressure 
exponent. 

3.2  APPLICATION  OF  SIMPLIFIED  "P-DOT"  THEORY 

Where  there  are  differences  between  closed  vessel  and  strand  data,  such 
differences  have  been  attributed  to  the  deductive  nature  of  the  closed  vessel 
experiment.  The  closed  vessel  experiment  measures  a variable  pressure  versus 
time  whereas  the  strand  experiment  measures  directly  the  distance  burned  versus 
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time  at  a constant  pressure.  From  the  pressure  versus  time,  ballistics  analysis 
yields  mass  generation  rate  versus  time  if  certain  assumptions  are  made.  If  the 
geometry  of  the  sample  is  known,  the  burning  rate  versus  time  may  then  be  calculated. 
Thus,  burning  rate  may  be  cross-plotted  versus  pressure.  However,  it  has  been 
suggested  that  the  assumptions  may  be  incorrect  and  that  the  geometry  of  the  sample 
may  not  be  known.  Unknown  factors  frequently  mentioned  include  heat  losses,  non- 
uniform  ignition  and  burning,  and  charge  break-up  (Ref.  18).  Low  pressure  uncertainties 
have  limited  the  reporting  of  data  to  pressures  above  5000  psi. 

A mechanistic  basis  for  a difference  between  burning  rates  measured  under 
steady-state  and  transient  pressure  conditions  was  developed  in  the  course  of 
studies  of  ignition  and  depressurization  transients  in  solid  rocket  motors  (Ref.  13). 
This  mechanism  was  the  thermal  wave  relaxation  time  in  response  to  a pressure  transient. 
A simplified  form  of  expression  was  as  follows  (Ref.  14): 


1 + 


2n< 


P 

P 


(1) 


r..  = instantaneous  burn  rate  at  pressure,  P 

r = steady-state  burn  rate  at  P 
n = steady-state  pressure  exponent  at  P 
k - propellant  thermal  diffusivity 
P = instantaneous  rate  of  change  of  pressure  at  P 


i 

I 

I 


More  detailed  analysis  resulted  in  the  development  of  computer  programs  such 
as  described  in  Ref.  (16).  Some  of  these  programs  have  been  reviewed  by  Kooker 
(Ref.  25)  in  the  context  of  gun  ballistics.  ' 
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An  order-of-magnltude  analysis  of  Eq.  (1)  reveals  that  P-dot  effects  should 
be  small  In  the  usual  closed  vessel  experiment.  Typical  values  of  the  parameters 
are  ad  follows: 

r = 2 In. /sec  @ 10,000  psi  (Fig.  1) 


n = 0.8  (Fig.  1) 

2 

< =0.0002  In.  /sec.  (typical) 


P = 1000  P (Ref.  20) 


The  Impact  of  P-dot  upon  burning  rate  falls  below  10%  at  pressures  above 
10  Kpsi , and  decreases  with  Increasing  pressures  and  burning  rates.  This  result 
Is  consistent  with  the  data  In  Fig.  1.  Although  the  values  of  P-dot  are  very 
high  when  compared  to  values  encountered  In  rocket  motors,  the  values  of  pres- 
sure and  r-squared  also  become  very  high  so  as  to  diminish  the  effect.  Since 
most  of  the  reported  closed  vessel  data  are  at  such  high  pressures,  it  would 
appear  that  agreement  between  the  closed  vessel  and  strand  data  should  be 
expected  assuming  that  P-dot  would  be  the  only  reason  for  any  difference.  This 
appears  to  be  borne  out  by  the  homogeneous  propellants. 

Results  of  P-dot  calculations  applied  to  the  nltramine  propellant  data  are 
Included  as  the  dashed  lines  In  Fig.  2.  It  was  possible  to  do  so  because  the 
P-dot  histories  for  these  tests  were  available  (Ref.  20).  Note  that  the  calcu- 
lated effect  has  the  tendency  to  reduce  the  pressure  exponent  In  the  high  slope 
region  below  5 Kpsi.  The  calculated  effect  becomes  small  at  high  pressure,  as 
discussed  previously.  However,  the  data  behave  quite  differently  from  the  cal- 
culated effect.  In  view  of  the  fact  that  experimental  uncertainties  and  P-dot 
effects  appear  to  be  adequately  resolved  In  the  case  of  homogeneous  propellants. 
It  follows  that  something  different  must  be  happening  In  the  case  of  nltramine 
propellants. 
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3.3  THE  ANALYTICAL  MODEL 

3.3.1  Conceptual  Basis 

| i 

The  analytical  model  for  the  transient  combustion  of  nitramine  propellants  ' \ 

is  based  upon  the  nitramine  surface  structure  mechanisms  contained  within  the  j 

steady-state  model  (Refs.  1,  7),  coupled  to  a thermal  wave  relaxation  model 

(Ref.  16)  of  a homogeneous  solid.  j 

! 

At  low  pressure,  and  following  ignition,  a planar  melt  layer  must  form  ! 

j ; 

and  develop.  This  development  will  be  an  endothermic  contribution  to  the  j 

energy  balance,  will  tend  to  oppose  the  ordinary  collapsing  of  the  thermal  wave  j : 

1 

during  a pressure  rise,  and  will  tend  to  delay  a break  point.  These  effects 
should  produce  a burning  rate  lag  at  low  pressure.  At  high  pressure,  and  where 
coarse  particles  are  present,  the  development  of  the  surface  craters  in  accordance 

1 

t 

with  break  points  will  increase  the  burning  surface  area  and  thereby  augment  the  j 

mass  burning  rate.  As  discussed  earlier,  the  closed  vessel  deduces  the  mass 
burning  rate  and  does  not  measure  directly  the  linear  burning  rate.  Therefore, 
if  the  burning  surface  is  underestimated,  the  linear  rate  as  computed  from  the 
closed  vessel  data  will  appear  to  lead. 

Examples  of  cratered  surfaces  appear  in  Fig.  9 of  Ref.  (1).  If  the  cratering 
mechanism  increases  both  the  linear  burning  rate  and  the  surface  area  of  the 
small  pellets  used  in  transient  devices,  then  the  mass  burning  rate  would  be 
proportional  to  the  square  of  this  process.  Thus,  the  exponent  break  point 
would  have  a more  serious  effect  in  pressurizing  chambers. 

3.3.2  Basic  Equations  for  a Homogeneous  Propellant 
Instantaneous  burning  rate  is  related  to  a reference  or  steady-state  value 

by  means  of  an  Arrhenius  equation  for  surface  decomposition: 
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r1  * r exp  [<-Eg/R0)(l/Tg  - 1/TS)  ] 


(2) 


[ I 


E$  * activation  energy  of  surface  decomposition 
Rq  * universal  gas  constant 
T$  = Instantaneous  surface  temperature 
Tg  = surface  temperature  corresponding  to  r 
r = steady-state  burning  rate  at  P,  known  from  data  or 
the  steady-state  model 


The  Instantaneous  surface  temperature  Is  calculated  from  the  time -dependent 
Fourier  equation: 


3 T 
3 t 


32T 


3X 


' "1 


Jl 

3X 


AQ 

-T^exP  (-ES/R0T) 


(3) 


T = temperature  distribution 
t = time 

x * distance  Into  the  solid 
Qss  * condensed  phase  heat  of  decomposition 
A = prefactor 

c = propellant  heat  capacity 

For  numerical  solution,  Eq.  (3)  Is  written  In  a suitable  finite  difference 
form  such  that  the  nominal  mesh  size  Is  selectable  depending  upon  the  thermal 
wave  thickness.  The  selection  affects  accuracy,  stability  and  computer  time. 
The  Initial  and  boundary  conditions  for  Eq.  (3),  using  Ref.  (1)  as  the  basis 
for  an  energy  balance,  are: 


10 


T (x,0)  - Tq 
£ <»•*>  ■ 0 


Qs-k  g(O.t)  • fext  ♦ f-1  Ks  - V * Oss]  e 


w = web  thickness 

Tq  = initial  bulk  temperature 

f . = external  heat  source,  for  example  as  would  be 
exz  Input  to  produce  Ignition 

k = propellant  thermal  conductivity 

Qs  = heat  flux  at  the  propellant  surface 

p = propellant  density 

£ = dimensionless  flame  height  corresponding  to  r at  P 

The  pressure  variations  are  computed  by  a mass  and  energy  balance  for  the 
chamber: 

(v-nmc)  ^ = (RTg  + "p)  Sf  + Rmc  "d^  “ p 31 


Tg  = gas  temperature 
mc  = mass  of  gas  stored  in  chamber 
R = gas  constant 
V = chamber  gas  volume 
n = covolume 


dmc 

dt 


dV 

’Ab  ri  and  dt  = Abr1 


= burning  surface  area 


For  a closed  vessel,  it  is  assumed  that  there  is  no  mass  exhaust. 
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dT 


dm. 


in 


_1  « (yT  - T ) -£■  - -*■ 
dt  'Y  F V dt  Cg 


(7) 


y ■ ratio  of  specific  heats 

Tp  ■ adiabatic  flame  temperature  (a  function  of  P) 

cg  ■ gas  heat  capacity  (a  function  of  P) 

* heat  loss  rate  to  wall 


1 

\ 


V\  ■ hc<Tg  - V + <“>,  <Tg4  - O <«> 

A^  * wall  area 

hc  * suitable  heat  transfer  coefficient 
* wall  temperature 

(eo)g  a product  of  gas  emlssivlty  and  Stefan  - Boltzmann  constant 

The  unsteady  heat  transfer  In  the  wall  Is  solved  in  analogous  fashion  to 
Eqs.  (3)  and  (4),  using  Eq.  (8)  to  describe  the  flux  at  the  boundary;  but  r^=0  and 
Q =0  for  the  wall. 

The  foregoing  equations,  plus  Eq.  (3)  written  for  the  wall,  constitute  a system 

• • 

for  the  unknowns  P,  V,  mc»  Tg,  Qw>  Tw»  Q$,  Tg  and  r^  as  functions  of  time.  These 
equations  are  solved  numerically  by  a forward  difference  method.  The  instantaneous 
burning  rate  will  differ  from  the  steady-state  value  at  a given  pressure  to  the  extent 
that  T$  differs  from  T . This  difference  can  arise  as  a result  of  the  changing  thermal 
profile  In  the  solid  In  response  to  the  transient  heating. 

3.3.3  Equations  Related  to  the  Surface  Structure  of  Nltramlne  Propellants 

Additional  equations  enter  Into  the  system  for  nltramlne  propellants  because 
of  the  relationship  of  the  surface  structure  to  the  burning  rate  (Refs.  1,  7)  and 
mass  flow  rate.  It  Is  assumed  that  the  surface  structure  will  develop  In  three 
stages  In  the  course  of  the  pressurization  transient: 

i 

i 

r 
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a)  formation  of  the  melt  layer  following  Ignition; 

b)  disappearance  of  the  melt  layer  when  the  break  point  criterion  Is 
reached; 

c)  formation  of  the  surface  craters  following  disappearance  of  the  melt 
layer. 

The  melt  layer  will  start  to  form  as  soon  at  = T^,  the  melting  point  of 
the  nitramine.  The  significance  of  the  melt  layer  thickness  is  that  it  positions 
QlM.  the  heat  of  fusion,  which  subtracts  from  the  exothermic  Q$$  in  the  energy 
bala  ''.e  equations.  It  is  expected  that  the  location  of  TM  will  move  into  the  solid 
with  time;  TM  is  551°K  for  HMX,  and  computed  steady-state  surface  temperatures  are 
circa  1000°K. 

The  second  stage  commences  when  the  critical  burning  rate  is  exceeded, 
r^>rm,  such  that  subsurface  melting  can  no  longer  propagate.  At  this  time  step, 

Qjj^l  is  set  to  zero  and  the  existing  melt  layer  thickness  is  allowed  to  decompose 
away  in  accordance  with  the  time  integral  of  r^  in  succeeding  time  steps.  This 
stage  is  completed  when  the  melt  layer  has  disappeared. 

The  formation  of  the  craters  impacts  both  r and  Ab  in  accordance  with  the 
steady-state  model.  However,  the  development  is  allowed  to  proceed  continuously 
with  time  (in  the  steady-state  model  it  appears  as  a discontinuity  with  pressure). 
The  changing  surface  area  in  the  course  of  crater  development  is  expressed  as 
follows: 


A^p  ■ "frozen",  or  cratered  surface  burn  area 


Ajj  = geometric,  or  input  burn  area 

(S  /S  ) = normalized  oxidizer  surface  exposed  for  a planar 

M surface,  from  Ref.  (7). 
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The  parameter  (S0X/SQ)p  changes  with  time  because  the  oxidizer  (nltramlne) 
penetration  depth,  h,  changes  with  time.  The  relationship  between  (s0X/S0)p 
and  h appears  In  Ref.  (1).  The  variation  of  the  penetration  depth  with  time 
may  be  expressed  as  follows: 


dh 

r 


(rox  - rf ) dt 


o*  ’ \w  exp  l-VW 


rf  = Af  exp  (-Ef/R0T$) 


(10) 

(11) 

(12) 


rox  = n^tram^ne  monopropellant  burning  rate 


rf  = binder  regression  rate 


Aqx  * kinetics  prefactor  for  the  nltramlne 


Eox  = act^vat^on  energy  for  the  nltramlne 


A^  = kinetics  prefactor  for  binder 
= activation  energy  for  binder 

A limit  Is  placed  upon  the  extent  of  penetration  In  accordance  with  the  steady- 
state  model.  This  limit  results  In  a maximum  (S0X/$0)p»  or  (s0X/s0)max*  From 
the  steady-state  model.  It  Is  known  that  the  area  change  from  (S  /S  )u  to 
(S0x/S0)f  Is  the  essential  manifestation  of  the  break  point  wherein  F jumps 
from  a value  Influenced  by  the  binder  to  a value  close  to  the  nltramlne  mono- 
propellant rate.  Therefore,  In  the  course  of  the  break  point  transition,  it 


may  be  assumed  that  r Is  proportional  to  the  change  in  (S0X/SQ).  This  assumption 


Is  made  In  order  to  avoid  jump  discontinuities  In  the  time-dependent  model,  and 
allow  smooth  (continuous)  transition.  Eq.  (2),  as  applied  during  this  third 
stage  of  surface  development,  becomes: 


i 

« s 


ZT  = exp  [(-Es/R)  (1/T$  - 1/T$) 


iL  , i.  , ^ox^F'^ox^M 
j _rM  . _(5oxV(S0xV 


(13) 


Fp*  post-break  steady-state  burn  rate 
rM«  pre-break  steady-state  burn  rate 
(5^)p  ■ that  value  associated  with  Fp 

14 


It  Is  anticipated  that  the  break  point  In  the  time-dependent  model  need  not 

occur  at  the  same  pressure  as  In  steady-state.  That  Is,  the  pressures  for 

r.>r  and  r>  r may  be  different.  Thus  It  is  necessary  to  include  values  of 
1mm 

r which  extrapolate  planar  surface  values  (r^)  to  higher  pressure  and  cratered 
surface  values  (Fp)  to  lower  pressure.  If  the  cratering  actually  occurs  at 
some  different  pressure  than  steady-state,  a proper  value  of  F may  still  be  used 
in  accordance  with  the  surface  existing  at  the  time. 

The  additional  unknowns  presented  by  nitramine  propellants  are  (S0X)M»  rQX, 
rf,  h,  (SQX)p  and  Abp.  Eqs.  (9-12)  and  two  equations  from  the  steady-state 
model  are  included  in  the  forward  difference  scheme  to  provide  the  relationships. 

Eq.  (13)  substitutes  for  Eq.  (2),  and  Abp  is  used  in  Eq.  (6).  The  changing  surface 
structure  is  an  added  means  by  which  the  instantaneous  burn  rate  may  differ  from 
the  equilibrium  value,  and  a difference  between  Ab  and  Abp  would  bear  upon  effective 
burning  rates  deduced  from  the  closed  vessel  experiment.  The  effective  burning 
rate  is  obtained  by  multiplication  of  Eqs.  (9)  and  (13). 

3.3.4  Computer  Programming 

A computer  program  providing  for  the  numerical  solution  of  the  equations  was 
written  for  operation  on  a UNIVAC  1108  computer.  The  program  has  been  successfully 
rendered  operational  and  applied  to  test  cases.  An  important  limitation  is  the 
run  time  which  can  be  encountered  at  very  high  burning  rates,  because  a rapid 
transient  demands  a smal  1 time  step  for  accuracy.  To  minimize  the  impact  of  this 
dependence,  a scheme  which  automatically  adjusts  the  time  step  and  the  solid 
propellant  grid  to  the  existing  transient  was  successfully  developed.  The  run 
time  for  a case  in  which  the  burning  rate  reached  20  in/sec  was  15  minutes.  A 
description  of  the  computer  program  is  presented  in  Appendix  B.  A source  deck 
and  listing  can  be  furnished  upon  request. 
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i 3.4  MODEL  RESULTS 

3.4.1  Homogeneous  Propellant  Test  Case 

The  homogeneous  propellant  test  case  selected  for  model  validation  is  "Test 
i Propellant  A"  as  specified  by  the  JANNAF  Workshop  on  Burning  Rate  Measurements, 

Memoranda  of  26  January  1977  and  5 April  1977.  This  case  was  selected  because  of 
the  extensive  data  provided  In  association  with  the  purposes  of  the  Workshop. 

Burning  rates  were  reported  by  six  different  laboratories  for  a given  set  of 
specifications  and  closed  vessel  pressure-time  data.  However,  the  identity  of 
the  propellant  was  not  disclosed  and  strand  burning  rates  were  not  given.  In  the 
absence  of  strand  data,  the  reported  burning  rates  were  used  as  baseline  equilibrium 
data  for  purposes  of  model  application.  The  data  were  reported  over  the  pressure 
range  5Kps1-40Kps1 . An  average  of  the  results  for  the  six  laboratories  was  used  at 
each  pressure,  and  an  extrapolation  was  performed  to  complete  the  input  down  to 
ambient  pressure.  Although  these  data  are  closed  vessel  data,  the  model  would 
still  be  useful  to  show  the  effect  of  the  closed  vessel  environment,  if  any.  If 
such  effect  is  small,  as  expected,  the  model  should  reproduce  these  data  as  well 
as  the  measured  pressure-time  curve. 

The  burning  rate  results  are  shown  In  Figure  3.  The  solid  line  represents  the 
actual  closed  vessel  data,  except  for  pressures  below  5Kps1  wherein  it  is  an 
extrapolation.  The  model  prediction  Is  represented  by  the  circles.  It  is  observed 
that  the  prediction  Is  very  close  to  the  line  between  5Kpsi  and  40Kpsi,  indicating 
a very  small  effect  of  the  transient  environment.  In  fact,  the  effect  of  the 
transient  environment  disappears  completely  for  this  test  case  above  13Kps1. 

An  Interesting  undulation  appears  In  the  predicted  results  between  4Kps1  and 
12Kps1.  It  Is  not  certain  whether  this  reflects  a real  mechanism,  a tradeoff 
between  the  transient  combustion  zone  heating  and  transpiration  cooling  during 
the  pressure  rise,  or  Is  simply  a result  of  the  forward  difference  numerical  method. 
Similar  undulations  have  appeared  In  computerized  data  reduction  for  closed  vessel 
burn  rate  (e.g..  Ref.  24). 
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Figure  3 


Predicted  Effect  of  Closed  Vessel  Environment  on  Burn  Rate 
of  a Homogeneous  Propellant 


The  most  interesting  behavior  appears  at  low  pressure.  The  very  high  initial 
burning  rates  are  largely  a consequence  of  the  assumed  Igniter.  The  model  does 
not  Include  a representation  of  the  Igniter  as  such,  but  uses  an  arbitrary 
Input  heat  flux  to  get  the  burning  under  way.  If  this  heat  flux  remains  operative 
during  burning,  and  Is  large  compared  to  the  combustion-induced  heat  flux,  the 
burning  rate  will  be  high.  The  Igniter  Is  turned  off  at  1000  psl,  and  the 
burning  rate  immediately  drops  because  the  transpiration  cooling  becomes  excessive 
for  the  remaining  heat  imposed.  The  burning  rate  climbs  again  when  the  transient 
heating  is  re-established.  Here,  then,  is  a second  type  of  burning  rate  undulation. 

A third  type  (not  shown  here)  results  from  pre-ignition  thermal  soak  of  the  propellant, 
or  dynamic  "overshoot"  due  to  a sudden  imposed  compression"  (25).  Thus,  there  are 
mechanistic  bases  for  low  pressure  burning  rate  undulations,  depending  upon  the 
nature  of  the  ignition  process  and  developing  combustion-induced  heating.  These 
effects,  although  real,  are  undesirable  for  purposes  of  closed  vessel  data 
acquisition  and  are  lumped  into  the  area  of  uncertainties.  As  a result,  burning 
rates  deduced  from  the  closed  vessel  generally  exclude  results  below  5 kpsi.  An 
accurate  quantification  of  the  low  pressure  behavior  would  require  a model  of  the 
igniter. 

Predicted  pressure-time  results  are  shown  together  with  data  in  Figure  4. 

Time  is  measured  from  ignition  time  rather  than  zero  time  to  remove  the  difference 
in  ignition  delay  and  afford  a better  visual  comparison.  Initial  pressurization 
is  overpredicted  because  of  the  effect  of  the  assumed  igniter.  This  also  causes 
the  predicted  web  to  advance  relative  to  the  actual  web  at  a given  time  and 
pressure.  The  curves  then  merge  and  the  prediction  becomes  very  good  at  the 
higher  pressures.  The  pressure  results  are  consistent  with  the  burning  rate  results. 


3.4.2  Nitramine  Propellant  Test  Case 


The  nitramine  propellant  test  case  selected  for  model  validation  is  the 
92-micron  HMX  propellant  reported  by  Lenschitz  (Ref.  20).  Much  of  the  required 


18 


78-6 


input  information  is  contained  within  that  paper,  including  actual  strand  burning 


rate  data  to  serve  as  the  equilibrium  baseline.  The  paper  also  reports  the 
pressurization  history  for  that  test.  Other  information  required  was  furnished 
by  Lenschitz. 

The  burning  rate  results  are  shown  In  Figure  5.  Plotted  are  the  strand  data 
(solid  line),  deduced  closed  vessel  data  (solid  circles),  and  the  predicted  actual 
(open  circles)  and  effective  (x)  burning  rates  in  the  closed  vessel.  As  expected, 
there  is  a significant  difference  between  strand  and  closed  vessel  burning  rates  in 
the  case  of  a nitramine  propellant,  and  the  model  serves  to  explain  why. 

First,  let  us  examine  the  region  below  2000  psi.  Here,  much  the  same  be- 
havior appears  as  was  discussed  for  the  homogeneous  propellant;  the  burning  rate 
is  first  dominated  by  the  behavior  of  the  assumed  Igniter,  and  drops  when  the 
igniter  is  turned  off.  Note,  however,  the  appearance  of  a second  burning  rate 
undulation  just  below  2000  psi.  This  is  a consequence  of  a fourth  mechanism  for 
undulation,  and  will  be  unique  to  nitramine  propellants  which  exhibit  a break 
point.  What  is  happening  here  Is  that  a break  point  is  being  traversed,  causing 
a greater  burning  rate  than  can  momentarily  be  sustained  by  the  combustion  heating 
then  in  existence.  The  rate  falls  until  sufficient  heating  develops  to  re-sustain 
it.  If  the  rate  fell  to  the  extent  that  the  surface  melt  layer  would  re-form, 
then  an  endotherm  would  reappear  to  drop  the  rate  further;  such  a process  might 
appear  as  a fifth  mechanism  for  undulation.  No  data  are  reported  below  2000  psi, 
but  the  model  results  are  not  expected  to  be  accurate  in  that  regime  because  the 
igniter  is  not  modeled. 

The  region  between  2000  psi  and  10  Kpsi  is  interesting  because  it  shows  how 
the  predicted  actual  burning  rates  in  the  closed  vessel  differ  from  equilibrium 
data.  Basically,  the  post-break-point  surface  structure  cratering  of  the  propellant 
is  developing  continuously  and  gradually  in  the  course  of  the  transient  in  the 
closed  vessel.  This  continuous  development  lags  the  step-wise  cratering  at  the 
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Figure  5.  Predicted  Effect  of  Closed  Vessel  Environment  on  Burn  Rate 
of  a Nltramlne  Propellant,  and  Comparison  With  Experimental 
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discrete  test-to-test  pressure  Intervals  In  the  strand  burner.  Thus,  the  severity 
of  the  break  point,  as  measured  by  the  exponent  In  this  transition  region.  Is 
lessened  by  the  closed  vessel  environment.  At  pressures  above  lOKpsI,  the  pre- 
dicted actual  burning  rates  In  the  closed  vessel  merge  with  the  equilibrium  strand 
data  as  the  surface  structure  becomes  fully  developed. 

Comparison  between  predicted  and  measured  burn  rates  In  the  closed  vessel  re- 
quires use  of  the  predicted  effective  burn  rate,  because  that  parameter  accounts 
for  the  effective  propellant  burn  surface  area  due  to  the  cratering.  The  predicted 
effective  burn  rate  does  not  depart  significantly  from  the  predicted  actual  burn  rate 
until  pressures  in  excess  of  3Kpsi.  Over  the  range  3Kpsi-20Kpsi , the  predictions 
are  generally  on  the  low  side,  but  agree  rather  well  with  the  data.  Thus,  It  is 
concluded  that  the  high  closed  vessel  burn  rates  at  high  pressure,  relative  to  the 
strand  data,  are  due  entirely  to  a surface  area  effect  and  not  a transient  burn 
rate  effect.  Lenschitz  had  rendered  a similar  conclusion  (Ref.  20),  but  thought 
that  the  propellant  grains  were  breaking  up.  Grain  breakup  is  not  required  to 
explain  these  results  if  the  cratered  surface  structure  is  recognized. 

The  predicted  pressurization  is  compared  with  data  in  Figure  6.  Here,  the 
data  were  reported  as  rate  of  change  of  pressure  versus  pressure  rather  than  as 
pressure  versus  time.  This  is  a more  severe  test  of  the  model  because  a number  of 
factors  enter  into  the  rate  of  pressure  change  at  a given  pressure.  The  prediction 
falls  below  the  data,  which  is  consistent  with  the  underprediction  of  the  effective 
burn  rates,  but  the  agreement  is  rather  good  considering  all  of  the  factors  involved. 

Based  upon  the  results  of  this  effort  and  the  steady-state  modeling  effort, 
it  would  be  expected  that  closed  vessel  burning  rates  should  agree  with  strand 
burning  rates  when  the  HMX  particle  size  Is  sufficiently  fine  that  there  is  no 
break  point  or  surface  cratering. 
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SECTION  4 

COMBUSTION  INSTABILITY  CHARACTERISTICS  OF  NITRAMINE  PROPELLANTS 

4.1  BACKGROUND  HYPOTHESIS 

A belief  that  HMX  may  have  favorable  effects  upon  combustion  Instability 
driving,  relative  to  AP,  Is  based  upon  the  qualitative  information  presented  in 
Table  1.  With  respect  to  some  factors  which  may  be  important,  HMX  appears  to 
be  more  similar  to  ammonium  nitrate  (AN)  and  potassium  perchlorate  (KP)  than  to 
AP.  Propellants  based  upon  AN  or  KP  have  not  been  reported  to  be  unstable 
(Ref.  26).  HMX  is  similar  to  AN  and  KP  In  that  there  is  a substantial  oxidizer 
melt  layer  and  particle  size  effects  upon  burn  rate  are  relatively  weak.1 
However,  HMX  surpasses  AP  as  an  energetic  monopropellant.  To  the  extent  that 
the  products  of  combustion  Include  particulates  which  may  afford  acoustic  damping, 
HMX  would  tend  to  produce  more  carbonaceous  residue  than  AN  (but  perhaps  not 
as  effective  as  spherical  KC1  droplets).  At  the  commencement  of  this  work, 
there  was  no  reported  systematic  study  of  the  effect  of  HMX  on  combustion  driving. 

4.2  PROPELLANT  SELECTION 

Standard  thermochemical  calculations  and  burning  rate  model  calculations 
were  performed  to  provide  selection  of  the  propellant  formulations  for  T-burner 
combustion  instability  experiments.  The  purpose  Is  to  determine  the  effect  of 
replacing  ammonium  perchlorate  (AP)  with  HMX  on  the  response  function  of  inert 
binder  (HTPB)  propellants,  and  of  HMX  addition  on  the  response  function  of 

nitrocellulose  (NC)  propellants.  In  making  these  comparisons.  It  was  desired 

* 

to  maintain  propellant  energy  and  burning  rate  reasonably  constant  within  each 
group.  Although  differences  can  be  corrected  in  the  course  of  interpreting  the 
response  function.  It  is  considered  that  excessive  differences  will  render  the 
comparison  less  meaningful  due  to  mechanistic  bias. 

T True  at  pressures  below  the  HMX  break  point.  Significant  changes  occur  above 
the  break  point,  which  would  alter  the  hypothesis.  Pressures  of  interest  In 
solid  rockets,  as  opposed  to  guns,  are  generally  below  these  break  points. 
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AP  and  HMX  produce  different  thermochemical  effects  In  HTPB  binder  because 
AP  Is  an  oxidizer  whereas  HMX  is  an  energetic  monopropellant.  HMX  propellants 
tend  to  be  low  burning  rate  propellants,  but  significant  burning  rate  control  Is 
available  In  AP  propellants  through  particle  size  adjustment.  Tradeoff  study 
produced  the  two  groups  of  selected  HTPB  propellants  listed  In  Table  2: 
a low  burning  rate  group  (XA)  and  a high  burning  rate  group  (XB).  It  is  difficult 
to  match  AP  propellant  burning  rates  with  HMX  alone,  so  the  burning  rate  of 
formulation  XA-3  is  lower  than  desired  for  Its  group  although  its  energy  level 
is  adequate.  The  two  mixed  oxidizer  propellants  are  more  closely  related  to  their 
AP  analogs  in  burning  rate  and  energy. 

It  was  found  that  the  addition  of  meaningful  quantities  of  HMX  to  an  inert- 
plasticized  NC  binder  produced  considerable  increases  in  propellant  energy.  Use 
of  a relatively  stable  energetic  plasticizer  (TMETN)  reduced  the  impact  of 
HMX  addition  on  propellant  energy,  and  was  considered  more  acceptable  for 
research  purposes  than  nitroglycerine.  The  resulting  propellants  comprise  the 
third  group  (XC)  listed  in  Table  2.  These  propellants  were  not  tested  in  the 
course  of  this  program  because  additional  testing  of  the  HTPB  propellants  was 
found  to  be  necessary  in  order  to  acquire  sufficient  meaningful  data. 
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4.3  EXPERIMENTAL  RESULTS 

Stability  tests  were  carried  out  in  2 1/2-inch  diameter  (I.D.)  T-burners 
coupled  directly  to  a surge  tank.  The  propellant  configuration  used  was  a 3/8- 
Inch  thick  disk.  Data  were  obtained  using  the  growth-decay  technique  (Ref.  27). 

For  each  propellant,  tests  were  conducted  at  nominal  frequencies  of  500,  900, 
and  1900  Hz  and  at  pressures  of  500  and  1000  Psi.  Several  tests  were  performed 
at  each  condition.  Many  of  the  tests  did  not  produce  oscillations,  or  the  data 
did  not  conform  with  T-burner  standards,  and  those  tests  are  excluded  from  the  results. 
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TABLE  2 

PROPELLANT  FORMULATIONS  FOR  T-BURNER  TESTS 


1 

XA-1 

XA-2 

XA-3 

XB-4 

XB-5 

XC-6 

XC-7 

i Wt-Z  AP(9tO 

- 

- 

- 

41.25 

31.25 

- 

- 

| * AP(50y) 

40.00 

20.00 

- 

- 

- 

- 

- 

AP(90p) 

- 

- 

- 

41.25 

31.25 

- 

- 

t * AP(200vj) 

| 

40.00 

20.00 

- 

- 

- 

- 

- 

\ HMX(4y) 

- 

20.00 

40.00 

- 

22.50 

- 

50.00 

1 HMX(15p) 

i 

- 

20.00 

40.00 

- 

- 

- 

- 

i 

i R-45/HTPB 

f 

20.00 

20.  OC 

20.00 

17.50 

15.00 

- 

- 

s 

| NC/TMETN 

- 

- 

- 

- 

- 

100.00 

50.00 

\ (T„/M) 1/2 

jf  . * 

10.34 

10.30 

10.32 

10.70 

10.77 

10.49 

10.86 

j (°K)1/2 

S r1000(in/sec) 

4 

0.261 

0.211 

0.161 

0.481 

0.431 

0.30 

0.33 

J 1 - Measured  Values 

I 

| 

| * 

1 

| * 

I 

l * 
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The  results  for  the  low  rate  and  high  rate  HTPB  propellants  are  given  In 
Figures  7 and  8,  respectively.  The  results  are  plotted  as  the  sum  of  the  measured 
growth  and  decay  coefficients  vs.  frequency.  The  propellants  were  all  relatively 
stable  In  that  the  measured  growth  coefficients  exceeded  the  decay  coefficients 
only  for  the  two  high  rate  propellants  at  the  1900  Hz,  500  psl  test  condition. 

Such  behavior  may  subject  the  absolute  results  to  criticism,  but  the  relative 
results  are  considered  valid. 

The  low  rate  XA-1  propellant  was  the  weaker  acoustic  driver  of  the  two 
propellants  which  did  not  contain  HMX.  The  low  rate  results  from  the  use  of 
coarser  particle  sizes.  Its  results  appear  to  be  relatively  independent  of 
pressure.  The  minimum  between  500  Hz  and  1800  Hz  is  considered  real  because  only 
one  test  out  of  seven  conducted  at  900  Hz  produced  measurable  acoustic  driving, 
whereas  most  of  the  excluded  data  at  the  other  frequencies  stemned  from  poor 
raw  data. 

The  two  low  rate  propellants  containing  HMX  exhibited  unusual  results.  At 
500  Hz,  the  pressure  traces  showed  an  immediate  growth  of  oscillations  following 
ignition  to  amplitudes  of  about  2%  of  the  mean  and  then  a slow  decay  over  the 
duration  of  the  burn.  At  1900  Hz,  the  tests  exhibited  no  oscillations  at  all. 

At  900  Hz,  the  500  psl  tests  were  like  those  at  500  Hz  and  the  1000  psi  tests 
were  like  those  at  1900  Hz. 

Post-test  examination  of  the  XA-3  tests  showed  the  T-burner  to  be  full  of  a 
layered,  carbonaceous  char.  In  the  1900  Hz  tests,  this  char  appeardd  to  have 
been  pulverized.  For  the  lower  HMX  concentration  XA-2  propellant,  the  char  was 
less  extensive,  merely  coating  the  wall  of  the  T-burners  rather  than  filling  the 
interior.  The  char  is  presumed  to  be  due  to  the  fuel -rich  nature  of  the  two  low 
rate  propellants  containing  HMX.  It  was  somewhat  surprising  because  prior 
experience  with  similar  propellants  did  not  encounter  char.  A review  of  the 
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) formulations  revealed  that  the  prior  propellants  were  much  more  heavily  plasticized 

igun  propellants,  whereas  the  propellants  formulated  here  were  more  consistent 
( with  the  structural  integrity  requirements  of  rockets.  Another  consequence  of 

I the  reduced  plasticizer,  or  the  char,  was  that  the  burning  rate  was  about  20% 

lower  than  expected.  The  char  is  indicative  of  non -equilibrium  thermochemistry 
in  that  there  should  be  no  free  carbon  at  a solids  loading  of  80%  HMX.  In  any 
event,  the  contribution  of  HMX  to  stability  could  not  be  Identified  between 
combustion  driving  or  mechanical  suppression  from  these  tests. 

The  high  rate  propellant  data  in  Fig.  8 Is  more  conclusive  as  to  the  role 
of  HMX.  These  propellants  did  not  produce  char.  Both  XB-4  and  XB-5  showed 
clear  pressure  and  frequency  effects,  increasing  In  driving  strength  with  de- 
creasing pre<>.  ' and  increasing  frequency.  The  XB-5  HMX-containing  propellant 
would  not  drive  oscillations  at  all  at  the  1000  psi  - 500  Hz  test  condition, 
and  its  data  at  500  psi  and  1000  psi  fall  below  the  respective  data  for  the 
all-AP  oxidizer  XB-4  propellant.  The  measured  damping  coefficients  at  each 
frequency  were  about  equal  for  the  two  propellants,  indicating  that  the  differences 
in  acoustic  driving  are  real.  More  data  are  required  for  clean-burning  propel- 
lants, and  testing  of  the  XC  group  is  included  in  future  plans. 
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SECTION  5 


EFFECT  OF  NITRAMINE  INGREDIENTS  ON  STEADY-STATE  BURNING  RATE 
(Data  Supplement  to  Annual  Progress  Report  AF0SR-TR-76-il36) 

5.1  PREFACE 

Strand  burning  rate  data  supplementing  the  data  reported  In  the  FY  1976 
Annual  Progress  Report  (Ref.  1)  were  furnished  by  AFATL  following  the  publica- 
tion of  that  Report.  These  data  are  presented  and  discussed  In  the  following 
subsections.  It  is  concluded  that  the  data  provide  further  confirmation  of  the 
results  and  predictions  of  the  steady-state  model  developed  in  the  course  of 
that  research. 

5.2  EFFECT  OF  HMX  PARTICLE  SIZE  IN  ACTIVE  BINDER  PROPELLANTS 

Figure  5 of  Ref.  (1)  showed  the  effect  of  HMX  particle  size  on  the  burning 
rate  of  nitroplasticized  polyurethane  (NPPU)  propellant.  Analogous  data  for  NC 
propellant  are  now  available  and  are  presented  in  Fig.  9.  The  same  basic  trends 
appear'  as  in  the  NPPU  propellant,  as  was  expected.  The  fine  propellant  exhibits 
a relatively  uniform  pressure  exponent  over  two  orders  of  magnitude  of  pressure. 

The  coarse  propellant  begins  to  depart  significantly  at  a pressure  above  lKpsI, 
with  an  upward  break  in  pressure  exponent  covering  the  range  1.3Kpsi  to  2.5Kds1. 

A second  upward  break  appears  to  begin  at  4Kpsi  and  continue  to  high  pressure. 

As  a Result,  the  burning  rates  of  the  coarse  propellant  have  a multiple  zig-zag 
appearance  and  are  higher  at  high  pressure.  The  burning  rates  of  the  NC  propellants 
are  higher  than  those  of  their  NPPU  analogs,  which  is  attributed  to  the  higher 
energy  and  burning  rate  of  the  NC  binder. 

A similar  trend  is  shown  by  the  mixed  HMX-TAGN  propellants  in  Fig.  10.  A 
portion  of  the  data  for  the  propellant  containing  the  fine  HMX  appeared  in  Fig.  7 
of  Ref.  (1).  No  data  were  available  for  the  propellant  containing  the  coarse 
HMX.  Results  for  both  propellants  are  now  complete.  The  propellant 
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containing  fine  HMX  exhibits  a relatively  uniform  shape  over  a wide  pressure 
range.  The  propellant  containing  coarse  HMX  exhibits  an  upward  break  In  pres- 
sure exponent  covering  the  range  lKpsi  to  3Kps1,  such  that  Its  burning  rates  are 
higher  at  high  pressures,  as  was  expected.  Comparing  Figs.  9 and  10,  the  effect 
of  the  TAGN  Is  to  raise  burning  rates,  more  so  at  low  pressure  than  at  high 
pressure,  such  that  the  overall  net  pressure  exponent  is  reduced.  This  effect 
has  been  attributed  to  the  faster  kinetics  of  TAGN  decomposition  and  combustion 
as  a monopropellant  (Ref.  1). 

Figs.  9 and  10  illustrate  the  effectiveness  of  fine  particle  size  as  an 
approach  to  achieve  desirable  slope  uniformity.  Fig.  9 confirms  that  it  is 
possible  to  do  so  with  HMX  alone  if  the  particle  size  is  sufficiently  fine. 

Fig.  10  confirms  that  TAGN  addition  is  helpful  provided  that  the  HMX  component 
does  not  exhibit  a slope  break. 


5.3  EFFECT  OF  TAGN  IN  ACTIVE  BINDER  PROPELLANTS 

Data  presented  in  Figs.  9 and  10  allow  the  completion  of  Fig.  7 of  Ref  (1). 
This  completed  figure  is  presented  as  Fig.  11.  The  effect  of  TAGN  substitution 
for  HMX,  in  increasing  amounts,  is  thereby  illustrated.  It  was  predicted  that 
the  data  for  the  29%  TAGN  propellant  would  interweave  with  the  data  for  the 
50%  TAGN  propellant  at  the  higher  pressures.  This  prediction  is  verified  by 
the  data.  A break  point  appears  for  the  50%  TAGN  propellant  at  lOKpsi.  This 
was  explained  by  the  burning  rates  reaching  such  a level  as  to  trigger  the  break 
in  the  HMX  component,  and  was  described  as  a disadvantage  of  excessive  TAGN 
addition.  In  view  of  that  explanation,  the  29%  TAGN  propellant  also  should 
exhibit  a break  point  when  the  burning  rates  reach  the  same  level.  This  is  not 
so  obvious  in  its  high  pressure  data,  perhaps  because  the  region  is  not  as  well 
characterized  with  data  points.  In  any  event,  it  is  confirmed  that  TAGN  addition 
raises  burning  rates,  and  that  optimum  amounts  are  helpful  as  an  approach  to 
achieve  desired  slope  uniformity. 
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There  has  been  a question  as  to  whether  or  not  propellants  containing  TAGN 
alone  would  exhibit  a break  point.  The  uppermost  data  in  Fig.  11  are  for  such 
a TAGN  propellant.  The  model  predicted  that  a break  point  exists  near  the  region 
where  considerable  data  scatter  appears,  and  the  scatter  was  Interpreted  to  be  a 
manifestation  of  the  break.  The  model  also  predicted  that  the  second  upward 
break  in  the  50%  TAGN  propellant  Is  due  to  the  TAGN.  The  question,  however, 
could  not  be  resolved  from  the  data.  Complete  data  are  now  available  for  the 
63%  TAGN,  NPPU  propellant  (cf.  Fig.  6 of  Ref  [1]).  These  data  are  shown  to- 
gether with  the  Fig.  11  data  in  Fig.  12.  A break  point  for  this  TAGN  propellant 
clearly  appears  at  lOKpsi.  Therefore,  it  can  be  concluded  that  TAGN  propellants 
can  exhibit  break  points.  Ironically,  this  break  point  can  be  mitigated  by 
lowering  the  burning  rates,  through  HMX  substitution,  which  would  yield  the  GAU-8 
propellant.  The  shape  of  the  NPPU  propellant  burning  rate  curve  lends  credence 
to  the  interpretation  of  a high  slope  (break)  region  in  the  case  of  the  NC 
propellant.  The  NC  propellant  has  higher  burning  rates  because  of  the  greater 
TAGN  concentration  and  the  more  energetic  binder. 

5.4  MATCHING  NITR AMINE  AND  ACTIVE  BINDER 

It  was  pointed  out  in  Ref.  (1)  that  matching  the  burning  rates  of  the  active 
binder  and  nitramlne  as  monopropellants  could  yield  a propellant  having  a uniform 
pressure  exponent.  The  rationale  was  that  a break  point  would  have  no  way  to 
manifest  itself,  and  so  would  not  appear  in  the  burning  rate  data.  In  support 
of  this  hypothesis,  data  and  results  of  model  calculations  were  presented  in 
Figs.  31-33  of  that  Report.  These  figures  showed  trends  as  the  nitramine 
burning  rate  approaches  and  falls  below  the  binder  burning  rate.  It  was  suggested 
that  a propellant  combining  fine  EDNA  and  a tailored  active  binder  having  the 
same  burning  rates  as  EDNA  would  be  a practical  test  of  the  approach. 
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Figure  12  Burning  Rates  of  TAGN  Active  Binder  Propellants 
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Unfortunately,  the  propellant  that  was  furnished  combined  coarse  EDNA  and 
the  baseline  NC  binder.  As  a result,  the  binder  burning  rates  are  lower  than 
the  EDNA  burning  rates  In  accordance  with  the  Figure  31,  and  the  break  point 
mechanism  would  be  triggered  at  a low  pressure.  With  this  propellant,  however. 

It  Is  possible  to  complete  the  EDNA  propellant  portion  of  Fig.  31  and  these  data 
are  shown  In  Fig.  13.  The  premise  that  the  post-break  propellant  burning  rates 
follow  the  nltramlne  monopropellant  burning  rates  Is  confirmed.  Without  data  at 
lower  pressures,  where  a break  would  be  expected.  It  does  not  appear  to  matter 
that  the  binder  burning  rates  are  low.  A transition  or  jump  upward  to  the  EDNA 
line  is  not  seen  over  this  pressure  range,  such  that  this  EDNA/NC  propellant 
appears  to  achieve  a desirable  result. 

An  illustration  of  binder  tailorability  to  Implement  the  foregoing  approach 
is  presented  in  Fig.  14.  The  lower  dashed  line  Is  for  nitrocellulose  (NC) 
containing  an  inert  plasticizer.  The  upper  dashed  line  is  for  nitrocellulose 
containing  a very  energetic  plasticizer.  The  line  fitting  AFATL  data  is  for 
nitrocellulose  con^'ning  an  energetic  (Intermediate)  plasticizer.  The  range  of 

lr 

burning  rates  is  observed  to  cover  a factor  of  3.  A larger  range  would  be  available 
by  adjusting  ingredient  proportions,  subject  to  energy  and  processing  limitations; 
a practical  range  of  5 is  indicated  by  existing  active  binders.  The  Intermediate 
binder  shown  in  Fig.  14  has  burning  rates  quite  close  to  EDNA  burning  rates,  and 
would  be  an  attractive  candidate.  Note,  however,  that  the  details  of  the  binder 
data  Indicate  some  break  point  mechanism  In  the  binder  which  might  be  reflected 
by  the  propellant  using  It.  It  Is  difficult  to  comment  upon  break  points  In  active 
binders  In  general  because  the  data  are  usually  reported  as  straight  lines. 
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Figure  13^  Comparison  of  Monopropellant  Burning  Rates  and  EDNA 
Propellant  Burning  Rates 
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1 SECTION  7 

| CONCLUSIONS 

A comparison  of  strand  and  closed  vessel  burning  rate  data  for  homogeneous 
propellants  reveals  that  effects  of  pressurization  are  minor,  and  that  the  two 
methods  should  give  consistent  results  at  pressures  following  the  ignition 
transient  in  the  closed  vessel.  A description  of  burning  rates  during  the 
ignition  transient  in  a closed  vessel  would  require  a model  of  the  igniter. 

The  burning  rates  of  nitramine  propellants  as  measured  in  a closed  vessel 
will  differ  from  strand  data  for  those  propellants  which  exhibit  a break  point. 

The  effect  of  the  transient  environment  in  a closed  vessel  is  to  create  a time 
lag  in  the  development  of  the  post-break-point  cratered  surface  structure,  which 
results  in  a less  severe  exponent  shift  than  shown  by  the  strand  data.  The 
reduced  severity  of  the  exponent  shift  will  have  a practical  benefit  in  transient 
devices.  However,  the  cratered  surface  structure  will  maintain  high  mass  burning 
rates  and  exponents  because  of  the  i ncreased  burn  surface  area.  Methods  to  avoid 
or  mitigate  the  cratering,  as  discussed  in  the  Annual  Progress  Report  for  FY  1976 
(Ref.  1),  are  still  recommended.  Additional  data  have  been  presented  here  which 
confirm  the  conclusions  of  that  Report. 

Several  mechanisms  have  been  identified  which  might  explain  burning  rate 
undulations  often  observed  from  closed  vessel  data  reduction,  in  particular  at 
low  pressure:  ignition  overshoot,  termination  of  the  igniter,  tradeoff  between 
developing  combustion  zone-induced  heating  and  increasing  surface  transpiration 
plus.  In  the  case  of  nitramine  propellants,  break-point  overshoot  and  melt  layer 
reformation. 

HMX -containing  composite  propellants  are  more  stable  than  comparable  AP 
propellants.  One  mechanism  may  be  related  to  mechanical  suppression  due  to  the 
generation  of  solid  material  combustion  products  in  certain  fuel -rich  formulations. 
However,  It  also  Is  true  that  HMX  reduces  the  combustion  driving  In  clean-burning 

i 

formulations  tested  thus  far.  Additional  study  is  recommended. 
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APPENDIX  A 
NOMENCLATURE 


n 

P 

P 

?ss 

% 

r 

rJ 

rF 


R 

Ro 

($ox/Vf 

^^F 

t 

T 


kinetics  prefactor  for  propellant  decomposition 
propellant  burn  surface  area 

nltramlne  propellant  burn  surface  area  following  cratered  surface  development 

kinetics  prefactor  for  binder  decomposition 

kinetics  prefactor  for  nltramlne  decomposition 

surface  area  of  closed  vessel  wall  exposed  to  heating 

propellant  heat  capacity 

gas  heat  capacity 

activation  energy  for  propellant  decomposition 
activation  energy  for  binder  decomposition 
activation  energy  for  nitramine  decomposition 
igniter  heat  flux 

depth  of  surface  cratering  In  nitramine  propellants 

heat  transfer  coefficient  for  wall  heatinq 

propellant  thermal  conductivity 

mass  of  gas  In  chamber 

burning  rate  pressure  exponent 

pressure 

rate  of  change  of  pressure 

heat  flux  at  propellant  surface 

propellant  heat  of  decomposition 

rate  of  heat  loss  to  closed  vessel  wall 

steady-state  burning  rate  at  P 

regression  rate  of  binder 

steady-state  post  break  point  burn  rate  at  P 

Instantaneous  propellant  burn  rate 

steady-state  pre-break  point  burn  rate  at  P 

regression  rate  of  nitramine 

gas  constant 

universal  gas  constant 

fraction  of  exposed  nltramlne  surface  on  a planar  melt  surface 
fraction  of  exposed  nitramine  surface  on  a cratered  surface 

value  of  associated  with  rr 

OX  r 

time 

temperature 
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w 

x 


e 

n 

Y 

K 

P 

<r 
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propellant  adiabatic  flame  temperature 

Instantaneous  gas  temperature 

Initial  (ambient)  temperature 

Instantaneous  propellant  surface  temperature 

steady-state  propellant  surface  temperature  associated  with  r 

closed  vessel  wall  temperature 

Instantaneous  gas  volume  In  the  closed  vessel 

propellant  web 

distance  from  propellant  surface  Into  propellant 

radiation  emlsslvlty  of  gas 
covolume 

ratio  of  specific  heats 
propellant  thermal  dlffuslvlty 
propellant  density 
Stefan-Bol tzmann  constant 

steady-state  flame  height  associated  with  r at  P (dimensionless) 


r 
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APPENDIX  B 


COMPUTER  PROGRAM  DESCRIPTION 


B-l  SUBROUTINES 

The  computer  program  file  is  identified  as  KT,  and  its  map  of  linked 
elements  for  execution  is  identified  as  KRTT.  The  component  subroutines 


are  listed  as  follows: 


Segment 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


Subroutine 


MAIN  PROGRAM 


SETUP 


DIFEQ 


TSWALL 


TSPROP 


PRINT 


A schematic  of  the  computer  program  elements  is  shown  in  Fig.  B-l. 

Segment  MAIN  is  the  master  control  element.  This  element  reads  input, 
computes  certain  parameters  of  an  input  nature,  writes  the  input,  and 
then  transfer  control  to  EVES.  Subroutine  EVES  is  a simultaneous 
differential  equation  solver  that  uses  the  Adams  predictor -corrector 
method.  EVES  sets  the  time  step,  predicts  ahead,  compares  the  results 
and  certifies  accepted  answers.  Subroutine  SETUP  initializes  all  of  the 
variables,  and  is  the  first  subroutine  called  by  EVES.  Thus  SETUP  provides 
the  initial  conditions  for  the  time-dependent  problem.  Subroutine  DIFEO 
computes  the  time  derivatives  of  gas  pressure,  gas  volume,  gas  weight, 
gas  temperature  and  distance  burned.  DIFEQ  contains  the  time-dependent 
conservation  equations  which  compute  these  differentials.  EVES  calls 
DIFEQ  after  the  problem  has  been  initialized.  DIFEQ  calls  TSPROP  to 
obtain  the  propellant  combustion  contributions  to  the  conservation 


B-l 


a)  MAX.  RUN  TIME  EXCEEDED 
OR  b)  MAX . PRESSURE  EXCEEDED 
OR  c)  WEB  BURN-OUT  ACHIEVED 


Figure  B-l.  Schematic  of  Computer  Program  Elements 
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equations,  and  TSWALL  to  obtain  the  heat  loss  to  the  wall  which  enters 
into  those  equations.  TSPROP  computes  the  time-dependent  temperature 
profile  in  the  solid  propellant,  the  burning  rate  which  is  a function 
of  surface  temperature,  and  includes  mechanisms  by  which  nitramine  behavior  can 
influence  the  burning  rate.  A finite  difference  approximation  of  the 
transient  heat  equation  is  used  to  compute  the  temperature  profile,  which 
develops  in  response  to  igniter-induced  heating  and  then  to  combustion 
heating.  TSPROP  may  control  the  time  step  when  high  burn  rate  leads  to 
steep  temperature  profiles.  TSWALL  computes  the  convective  and  radiation  heating 
of  the  wall,  and  the  temperature  profile  in  the  wall,  in  determining  the 
heat  loss.  TSPROP  and  TSWALL  return  answers  to  DIFEQ , and  DIFEQ  returns 
answers  to  EVES.  EVES  then  calls  PRINT  to  update  the  accepted  answers  and  to 
write  them  at  pre-selected  time  intervals.  Subroutine  FIND  is  an  auxiliary 
subroutine  that  is  used  to  interpolate. 

The  computations  involving  TSPROP,  which  is  the  essence  of  the 
combustion  model,  warrant  a more  detailed  discussion  which  is  provided  in 
Section  2 of  this  Appendix.  A source  deck  and  listing  of  the  computer 
program  can  be  furnished  upon  request;  it  is  liberally  embellished  with 
comment  cards  to  assist  the  user.  A manual  for  inputs  is  provided  as 
Section  3 of  this  Appendix.  Output  is  discussed  in  Section  4. 

B-2  THE  COMBUSTION  MODEL  CALCULATIONS  (TSPROP) 

The  following  sequential  discussion  describes  the  computations  performed 
in  TSPROP. 

1.  The  time  step  is  first  determined  by  EVES. 

2.  Tabular  strand  burning  rates  are  initialized  for  the  beginning  of 
this  time  step. 


78-6 


f 


I 


3.  The  propellant  web  calculational  mesh  and  counter  are  updated 
for  the  web  which  was  burned  away  in  the  previous  time  step. 

4.  A further  adjustment  to  the  mesh  is  made  depending  upon  whether 
the  reference  mesh  is  greater  than  1/10  the  thermal  wave  thickness 
of  the  previous  time  step,  or  less  than  1/25  of  that  thickness. 

Too  coarse  a mesh  impairs  accuracy,  too  fine  a mesh  will  increase 
computer  time  by  requiring  a smaller  time  step.  If  the  mesh  is 
adjusted,  the  temperature  profile  is  redefined  by  logarithmic 
interpolation  to  conform  to  the  new  mesh.  This  technique  of 
adjusting  the  mesh  has  provided  considerable  savings  in  computer  time. 

5.  The  time  step  required  by  the  mesh  is  compared  to  the  time  step 
required  by  EVES,  and  the  smaller  value  controls. 

6.  The  igniter -induced  heat  flux  is  gradually  diminished  to  zero  when 
sufficient  burning  has  been  achieved.  The  method  of  representing 
the  igniter  is  arbitrary.  There  is  no  model  of  the  igniter  as 
such,  only  an  input  heat  flux  to  start. 

7.  The  heat  flux  to  the  propellant  is  computed,  consisting  of  igniter 
and  combustion  contributions.  The  actual  heat  flux  is  referenced 
from  that  which  would  exist  under  equilibrium  burning  at  the  pres- 
sure of  the  present  time  step,  and  therefore  assumes  a quasi-steady 
gas  combustion  zone.  The  appropriate  tabular  strand  burning  rate  is 
updated  to  the  present  time  step  for  this  purpose.  If  the  propellant 
is  a nitramine  propellant,  and  the  actual  burning  rate  of  the 
previous  time  step  is  below  the  breakpoint  burning  rate,  then  the 
exothermic  surface  heat  release  term  is, modified  by  the  heat  of 
fusion  of  the  nitramine. 
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8.  Using  a finite-difference  expression,  the  resulting  heat  flux 
is  converted  into  an  artificial  surface  temperature  that  is 
used  to  drive  the  thermal  profile  in  the  depth  of  the  solid. 

9.  The  finite  difference  form  of  the  transient  heat  equation  is  used 
to  compute  the  temperature  at  each  mesh  location  in  the  depth  of 
the  solid,  including  the  true  surface  temperature.  The  rate  of 
change  of  temperature  is  computed  and  multiplied  by  the  time  step. 

10.  The  actual  burning  rate  is  computed  from  the  surface  temperature 
using  an  Arrhenius  expression. 

11.  Propellant  density  is  known  from  inputs.  The  present  burn  area 
is  known  from  subroutine  DIFEQ.  Knowing  the  actual  burn  rate 
the  propellant  weight  flow  rate  is  then  computed. 

12.  If  the  propellant  is  homogeneous,  this  burn  rate  and  weight 
flow  rate  are  returned  to  DIFEQ  at  this  point.  If  the  propellant 
is  a nitramine  propellant,  additional  computations  are  performed. 

13.  A test  is  made  comparing  the  actual  burning  rate  with  the  break -point 
burning  rate.  If  the  actual  burning  rate  is  less,  the  melt  layer 
thickness  of  the  nitramine  is  determined  by  the  location  of  the 
melt  temperature  in  the  thermal  profile.  If  the  actual  burning 

rate  is  now  less  but  was  at  one  time  greater  than  the  break  point 
rate,  such  that  the  surface  is  now  cratered,  then  an  additional 
computation  is  performed  to  begin  or  continue  to  knit  the  surface 
back  to  a planar  melt.  This  latter  condition  can  result  from  de- 
pressurization or  reductions  in  heat  flux. 

14.  After  computing  the  melt  layer  thickness,  the  program  returns  to 
DIFEQ  if  the  surface  is  now  or  is  once  again,  a planar  melt.  If  the 
surface  has  not  yet  knitted  completely,  then  adjustments  are  made 
according  to  the  extant  of  cratering  as  discussed  below. 
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| 15.  If  the  actual  burning  rate  is  greater  than  the  break -point 

« 

| burning  rate  in  step  13,  then  the  melt  layer  thickness  is 

reduced  by  the  product  of  the  burning  rate  and  the  time  step. 

If  the  melt  layer  thickness  is  still  greater  than  zero,  the  program 
returns  to  DIFEQ . If  the  melt  layer  has  disappeared  by  this 
mechanism,  then  the  extent  of  cratering  is  computed  according 
to  the  difference  between  the  nitr amine  and  binder  burning  rates. 
These  ingredient  burning  rates  are  calculated  from  Arrhenius 
expressions  for  the  present  surface  temperature.  If  the  binder 
is  an  active  binder,  its  rate  is  determined  from  the  updated 
tabular  input. 

'•  16.  If  the  surface  is  cratered,  either  as  a result  of  step  15  or 

step  14,  the  actual  burning  rate  is  upgraded  in  accordance  with 
the  mechanism  contained  in  the  steady-state  model.  The 
parameters  which  produce  this  upgrade  are  computed  from  the 
known  extent  of  cratering  and  tabular  burn  rates  as  presently 
updated.  Associated  with  this  upgraded  burn  rate  is  an  upgraded 
burn  area,  also  computed  from  the  cratering.  Thus  an  upgraded 
weight  flow  rate  is  computed,  and  these  upgraded  values  are 
returned  to  DIFEQ  instead  of  the  values  from  Step  12. 

B-3  CARD  INPUT  SEQUENCE  AND  STANDARD  VALUES 

The  following  describes  input  cards  in  sequence,  with  appropriate 
comments  and  numerical  values  which  have  been  used. 
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JOB  CARD 

SEQUENCE  OP  CONTROL  CARDS 
EXECUTION  CARD 


CARD 

COLUMN 

ENTRY 

SYMBOL 

UNITS 

1 

2 

Type  of  Propellant 

ICASE 

- 

Comment!  ICASE  = 1 for  homoqeneous  propellants 

ICASE  = 2 for  nitramine  propellants 

2 

3 

Title  of  Case  or  Group 

ITIT 

3 

2 

Initial  Chamber  Gas  Volume 

VINIT 

in  _ 

12 

Initial  Pressure 

PINIT 

lb/in 

22 

Initial  Temperature 

TINIT 

°K  2 

32 

Igniter  Heat  Flux 

FEXT 

cal/cm  -sec 

Comment:  Values  of  100-200  have 

been  used. 

42 

Reference  Calculational  Grid 

size  of  Propellant 

DELXO 

in. 

Comment:  A value  of  0.0002 

has  been  used. 

52 

Calculational  Grid  Size  for  Wall 

DELR 

in. 

Comment:  Divide  wall  thickness 

by  10 

3 

62 

Chamber  Diameter 

DIA 

in. 

4 

2 

Surface  Area  of  Chamber  Wall 

AWALL 

% 

cm  /sec. 

12 

Thermal  Diffusivity  of  Wall 

ALPHAN 

22 

Flame  Temperature  of  Propellant 

at  1000  psi 

TFP 

°K 

Comment:  Divide  isochoric  flame 

temperature  by  GAMMA 

32 

Gas  Molecular  Weight 

MW 

gm/mol 

42 

Ratio  of  specific  heats 

GAMMA 

- 

52 

Specific  heat  of  gas  at  10,000 

cal 

psi 

CGP 

gm-°K 

62 

Specific  heat  of  solid  propellant 

CSP 

cal 

5 

2 

Co volume 

COVOL 

gm-°K 
in  /lb. 

12 

Weight  Fraction  of  Nitramine 

ALFA 

- 

22 

Mean  Particle  Size  of  Nitramine 

DAVE 

Microns 

32 

Maximum  Surface  Penetration 
Associated  with  Post-Break -Point 
Surface  Cratering 

HDMAX 

Comment:  Use  a value  of  3.  for  sizes  greater  than 

100  microns 

Use  a value  of  1.  for  sizes  less  than 

10  microns 

Use  a value  of  2.  for  intermediate  sizes 

2 

42 

Nitramine  density 

RHOX 

gm/crn^ 

52 

Binder  density 

RHOF 

_gm/cm 

62 

Propellant  Thermal  Diffusivity 

ALPHAP 

in  /sec. 
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CARD 

6 


7 


COLUMN 


2 


12 


22 

32 


42 


52 


62 


2 


12 

22 

32 

42 


ENTRY  SYMBOL  UNITS  . 

Reference  Propellant  Surface  Temperature 

Temperature  TSREF  °K  I 

Comment:  For  homogeneous  propellants , select  j 

a value  consistent  with  the  burning  rate 

at  1000  psi  and  the  kinetics  expression  j 

for  surface  decomposition;  an  example  is  i 

523°K  for  0.35  in/sec.  and  an  activation  j 

energy  of  50000  cal/mol.  Surface  temperatures  • < 

for  homogeneous  propellants  are  in  this  , 

vicinity.  For  nitramine  propellants,  do 
the  same  but  with  a pre-break-point  burning 

rate;  for  example,  926°K  for  0.13  in/sec  and  * . 

50000  cal/mol.  Nitramine  surface  temperatures 

are  circa  900°K,  and  may  be  acquired  from  the  j 

steady-state  model.  j 

RREF  in/sec . 

, 

TMELT  °K 


Reference  burning  rate  associated  with 
TSREF 

Comment:  See  comment  for  TSREF 

Nitramine  Melting  Point 

Comment:  The  value  for  HMX  is  551. 

The  value  for  RDX  is  467. 


Activation  Energy  for  Nitramine  or  Homo- 
geneous Propellant  Decomposition  ESP 

Comment:  A value  of  50000.  has  been 
used  for  each. 

Activation  Energy  for  Binder 
Decomposition  EF 

Comment:  Use  a value  of  16900.  for 
HTPB.  For  active  binders, 
input  0. 

Prefactor  for  Nitramine  Decomposition  ASP 

Comment:  A value  of  0.5000  E + 10  has 
been  used  for  HMX  and  RDX. 

Prefactor  for  binder  decomposition  AF 


Comment:  A value  of  299  has  been  used  for 
HTPB.  For  active  binders,  input  0. 

Surface  Heat  Release  QSURP 

Comment:  A value  of  -225.  has  been 
used  to  represent  HMX  and  energetic 
homogeneous  propellants. 

Nitramine  Heat  of  Fusion  OLM 


Comment:  A value  of  132.5  is  used 
for  HMX  and  RDX. 

Nitramine  Thermal  Diffusivity 
Comments:  A value  of  0.0011  has 
been  used. 

Binder  Thermal  Diffusivity 
Comment:  A value  of  0.0011  has 
been  used. 

Nitramine  Thermal  Conductivity 
Comment:  A value  of  0.00049  has  been 
used. 


XKAPO 

XKAPF 

XLAMO 


cal/mol 


cal/mol 


cm  -sec 


gin 


cm  -sec. 


cal/gm. 


cal/gm 


cm  /sec . 


cm2/sec . 


cal 

cm-sec-  K 
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COLUMN 


ENTRY 


7 Cont’d  52  Binder  Thermal  Conductivity 

Comment:  A value  of  0.00044 
has  been  used. 

62  Wall  Thermal  Conductivity 

3 2 Number  of  Points  in  the  Table  of 

Burn  Area  vs  Web 

9 2 Web  Burned 

Comment:  The  first  table  entry 
should  be  zero 
12  Burn  area 

Comment:  This  is  the  total  burn 
area  at  zero  web  burned . 

REPEAT  FOR  NWEB  CARDS  TO  DEFINE  AREA  VS  WEB. 

NEXT  2 '’umber  of  Points  in  the  Table  of 

Pre-Break-Point  Burning  Rate 

vs.  Pressure 

NEXT  2 Pressure 

12  Pre-Break  Point  Strand  Burn  Rate 

Comment:  Cover  the  full  pressure 
range  with  the  series  of  cards, 
extrapolating  to  PINIT  and  to  the 
maximum  pressure  expected.  The 
series  may  be  omitted  for  homogeneous 
propellants,  whence  the  NPl  card 
should  be  omitted  also. 

REPEAT  FOR  NPl  CARDS  TO  DEFINE  PRE-BREAK-POINT 

BURNING  RATE  VS  PRESSURE 

NEXT  2 Number  of  Points  in  the  Table  of 

Post-Break-Point  Burning  Rate  vs 
Pressure 

NEXT  2 Pressure 

12  Post-Break-Point  Strand  Burn  Rate 

Comment:  Cover  the  full  pressure 
range  with  this  series  of  cards, 
extrapolating  to  PINIT  and  to  the 
maximum  pressure  expected.  This 
series  may  be  omitted  for  homogeneous 
propellants,  whence  the  NP2  card 
should  be  omitted  also. 

REPEAT  FOR  NP2  CARDS  TO  DEFINE  POST -BREAK-POINT 

BURNING  RATE  VS  PRESSURE 

NEXT  2 Number  of  Points  in  the  Table  of 

Active  Binder  Burning  Rate  vs 
Pressure 


SYMBOL 


XLAMF 


UNITS 


cm-sec-  K 


cm-sec-  K 


WEB1X 


lb/in. 2 
in/sec . 


lb/in. 

in/sec 


CARO 


COLUMN 

2 

12 
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ENTRY  SYMBOL 

Pressure  PBX 

Active  Binder  Strand  Burn  Rate  RBX 


UNITS 

lb/in2 
in/sec . 


Comment:  Cover  the  full  pressure  range  with  this 
series  of  cards,  extrapolating  to  PINIT  and  to 
the  maximum  pressure  expected.  This  series  may 
be  omitted  for  inert-binder  propellants,  but  is 
a must  for  homogeneous  propellants. 

REPEAT  FOR  NP3  CARDS  TO  DEFINE  ACTIVE  BINDER  BURNING  RATE  VS  PRESSURE 


LAST  CARD 
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B-4  OUTPUT 

The  initial  portion  of  the  output  displays  the  input  as  groups  labeled 
computation  parameters , gas  parameters,  propellant  parameters  and  the  tabular 
input  parameters.  The  answers  are  displayed  in  columns  consisting  of  time, 
web  burned,  pressure,  burn  rate,  propellant  surface  temperature,  melt  layer 
thickness,  crater  penetration  depth,  fraction  of  exposed  nitramine  surface  area, 
burn  area  and  effective  burn  rate.  Except  for  burn  area,  the  last  five  parameters 
are  associated  with  nitramine  propellants  only.  The  effective  burn  rate  is  the 
actual  burn  rate  multiplied  by  the  ratio  of  the  cratered  surface  burn  area  to 
the  ideal  geometric  burn  area,  and  would  be  the  value  of  burn  rate  deduced  from 
a closed  vessel  test  in  the  absence  of  knowledge  of  the  actual  (cratered)  burn 
area.  The  burn  area  which  is  output  is  the  geometric  (uncratered)  burn  area. 

Answers  cure  displayed  in  time  increments  of  0.1  msec. 

Computer  run  time  may  be  controlled  by  a time  stop  or  max  time  in 
accordance  with  particular  facility  usage,  and  is  usually  included  in  a control  card.. 
The  run  will  automatically  stop  when  the  web  burns  out  or  the  pressure  of  40000  psi 
is  exceeded.  This  maximum  pressure  appears  in  subroutine  PRINT,  and  could  be 
changed  if  desired.  Full  runs  have  ranged  from  5 min. -15  min.  CPU  on  a 
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